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Abstract. In this paper, through an exhaustive analysis within the Migdal-Eliashberg theory, we show the 
incompatibility of experimental data of RbsCeo with the basic assumptions of the standard theory of su- 
perconductivity. For different models of the electron-phonon spectral function a 2 F(fi) we solve numerically 
the Eliashberg equations to find which values of the electron-phonon coupling A, of the logarithmic phonon 
frequency Q\ a and of the Coulomb pseudopotential n* reproduce the experimental data of RbaCeo. We 
find that the solutions are essentially independent of the particular shape of a 2 F(fl) and that, to explain 
the experimental data of RbsCeo, one has to resort to extremely large couplings: A = 3.0 ± 0.8. This 
results differs from the usual partial analyses reported up to now and we claim that this value exceeds 
the maximum allowed A compatible with the crystal lattice stability. Moreover, we show quantitatively 
that the obtained values of A and J?i n strongly violate Migdal's theorem and consequently are incompat- 
ible with the Migdal-Eliashberg theory. One has therefore to consider the generalization of the theory of 
superconductivity in the nonadiabatic regime to account for the experimental properties of fullerides. 

PACS. 74.70.Wz Fullerenes and related materials - 74.20.-z Theories and models of superconducting state 
- 63.20.Kr Phonon-electron and phonon-phonon interaction 



1 Introduction 



In contrast to the cuprates, the other family of high-T c 
superconductors, the fullerene compounds, shows a quite 
more conventional phenomenology: their normal state 
properties are Fermi liquid-like, no stripe formation or 
signal of pseudogap appear above T c , they are s-wave su- 
perconductors with a sizeable carbon isotope effect |jj . It 
is certainly due to their apparently ordinary phenomenol- 
ogy that superconductivity in Ceo materials is now of- 
ten assumed to be correctly described by the conventional 
Migdal-Eliashberg (ME) theory of the electron-phonon 
driven superconductivity @||. The relatively high criti- 
cal temperatures of the A^Cqq fullerene compounds (up 
to T™ ax = 40 K for Cs 3 C 6 o under pressure 0]) are there- 
fore generally thought to be due to an optimized electron- 
phonon interaction achieved, according to the different 
proposed explanations, by a large electron coupling to 
the alkali phonons ||, to the Cgo rotational modes ||, or 
by the so-called curvature argument for the intramolec- 
ular couplings According to this latter theory, the 
electron-phonon coupling steadily increases with the cur- 
vature of the fulleride molecule so that compounds with 
smaller molecular radius (C36, C28) are expected to have 



critical temperatures even higher than those of the C60- 
based materials || . 

The above theories disregard however several aspects 
of the phenomenology of the fullerene compounds which 
do not fit into the standard ME scenario. In fact, like 
the high-T c copper-oxides, the fullerene compounds have 
extremely low charge carrier densities |lfj, have signifi- 
cant electron correlation p|,pT[, and are close to a metal- 
insulator transition [ [l2|JT^ | showing a strong dependence of 
T c upon doping and disorder Jl5| . From a ME point of 
view, these features tend to degrade the superconducting 
state so that it appears difficult to understand why the 
fullerene compounds should represent the best optimized 
ME materials. 

In this situation, it should be therefore of primary im- 
portance to assess to which extend the experimental data 
can be explained by the ME theory. This issue has been 
even more highlighted by the recent discovery of super- 
conductivity with T c = 52 K in FET hole doped C 6 o [@ ■ 
Until recently, however, the spread of experimental results 
has prevented a detailed and definitive analysis. Neverthe- 
less, in the last years, the experimental uncertainty for the 
alkali fullerene compound RbjCgo, which shows the high- 
est critical temperature at room pressure T c — 30 K within 
the A3C60 family, has been considerably narrowed. 
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In fact, resistive measurements on 99 % enriched 13 C 
single-crystals have recently established the carbon iso- 
tope coefficient ac with the best up-to-date accuracy, ac = 
0.21±0.012 17|, resolving therefore a long standing uncer- 
tainty on the value of ac in fullerides mainly due to partial 
isotope substitution and/or magnetic measurements on 
powder samples with broad transitions often wider than 
the isotope shift itself In addition, crossed tunnel- 

ing and optical transmission studies on the same sample 
have provided accurate measurements of the zero temper- 
ature energy gap 2 A leading to 2A/T C = 4.2 ±0.2 Q. 
This result is more accurate than previous spectroscopic 
studies and agrees with NMR and photoemission mea- 
surements [fl9|j20f| , moreover it has been obtained by both 
bulk- (optical transmission) and surface- (tunneling) sen- 
sitive measurements which have given often contrasting 
results. 

Although we shall briefly consider also other estima- 
tions existing in literature, in this paper we mainly focus 
on the above reported experimental data. The main point 
in fact is that these data are supplemented by quite small 
error bars, which permit a more rigorous analysis, We shall 
show however that different determinations of the data, in 
particular of the superconducting gap, would even more 
substain the results presented in this paper. 

The experimental values above discussed seem, at first 
sight, perfectly compatible with the ME theory. For ex- 
ample, T c = 30 K is not far from T c = 23.2 K of Nb 3 Ge 
and 2A/T C = 4.2 ± 0.2 is very close to 4.18 and 4.25 of 
VaGa and NbaAl(2), respectively Q. Finally, the carbon 
isotope coefficient ac = 0.21 is similar to the isotope ef- 
fects of the elemental ME superconductors Os (a = 0.2) 
and Mo (a = 0.33) p^| . The apparently ordinary values of 
T c , ac, and 2A/T C , independently considered, could there- 
fore be used as arguments in favour of the validity of the 
ME theory for RbsCgo- However, each of this values has 
little meaning if taken individually. For example, in Ref. 
|Jj the experimental data T c = 30 K and 2A/T C = 4.2, 
but not ac = 0.21, have been fitted by setting A = 1.16, 
fj,* = 0.1, and Q Xn = 302 K, while in Ref. Jl| T c = 30 
K and a c = 0.21, but not 2A/T C = 4.2, have been fit- 
ted by A = 0.9, n* = 0.22, and Q Xn = 1360 K. Although 
the discrepancies in the values of A and fi* are somehow 
acceptable, the two values of f2\ n differ by a factor of five. 

A previous partial analysis, based only on the values 
of the critical temperature T c — 30 K and of the iso- 
tope coefficient ac = 0.21, has pointed out an intrisic 
inconsistency of the ME framework with respect of the 
adiabatic assumption In this paper we extend this 
study by taking into account also the estimation of the 
superconducting gap 2A/T C = 4.2 and by considering re- 
alistic electron-phonon interactions related both to inter- 
and intra- molecular modes. We show that the experi- 
mental data of RbsCgo, when analyzed in the context 
of the fullerene compounds, T c = 30 K, ac = 0.21 and 
2A/T C = 4.2 are actually incompatible with the stan- 
dard ME framework. We accomplish this task in Sec. || 
by considering different models of the electron-phonon in- 
teraction and by numerically solving the ME equation in 



order to reproduce the experimental data. We find in Sec. 
|] that the resulting values of the electron-phonon interac- 
tion are always far too large to avoid lattice instabilities 
and to neglect the electron-phonon vertex corrections be- 
yond Migdal's limit, as required by the ME formulation 
[PP]. Finally, in Sec. || we propose that the main origin 
of the failure of the ME framework lies in the breakdown 
of the adiabatic hypothesis (Migdal's theorem) which, in 
doped fullerenes, is naturally driven by the small value of 
the Fermi energy. 



2 The Migdal-Eliashberg equations 

As pointed out in the introduction, the superconducting 
state of the fullerene compounds is quite often regarded in 
terms of the ME theory. Since the alkali doped fullerenes 
are three dimensional s-wave superconductors, the super- 
conducting properties are therefore thought to be cor- 
rectly described by the standard ME equations p,pi|: 
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In the above equations Z(iuj n ) = 1 — U(iu) n )/iuj n , where 
S(iuj n ) is the normal state electronic self-energy, A{iu>) is 
the Matsubara gap-function and uj n and ui m are Matsub- 
ara fermionic frequencies. 

In the above equations, a 2 F[f!) is the electron-phonon 
spectral function (also known as the Eliashberg function) 
which defines the electron-phonon coupling constant A and 
the logarithmic frequency J?i n through the following rela- 
tions [Ell: 



A 



2 / £a* (fl ), 



(3) 



\nf2 ln = jJ^-\n(n)a 2 F(n). (4) 

The dimensionless parameter /1 represents the effective 
Coulomb repulsion probed by the Cooper pair at the en- 
ergy scale lo c which is much larger than the phonon energy 
scale. It is clear that the value of the Coulomb parameter 
fi depends on the specific value of oj c which is somewhat 
arbitrary j24|. In fact a more sound quantity is the pseu- 
dopotential fi* defined as: 



1 + /l In (W c /i7max) 



(5) 



where J7 max is the maximum phonon frequency. The phys- 
ical properties of the superconducting state depend on \x* 
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Fig. 1. Plot of 2A/T C (solid line) and 2A /T C (dashed line) as 
function of A for an Einstein phonon spectrum and = 0. In 
inset, a zoom of the region around 2A/T C = 4.2. 



rather than \i or ui c [f24| . In the numerical solution of Eqs. 
(|l|)-(||) we have set oj c = 5J? max . Larger values of uj c do 
not modify the results and we have checked that different 
choises of uj c , although providing different values of \i as 
expected, lead to the same \i* via Eq. (g). 

The quantities we want to extract from the ME equa- 
tions (Q) and are the critical temperature T c , the zero 
temperature superconducting gap A and the carbon iso- 
tope coefficient ac- To achieve this, we must first choose 
the input quantities a 2 F(fl) and ^ (or /i*). Of course, the 
electron-phonon coupling A and the logarithmic phonon 
frequency fi\ n follows directly from a given Eliashberg 
function a 2 F(£2) via Eqs. (§)-(§)• Our aim is to find for 
which values of A, Sl\ n and /i* the ME equations have 
T c = 30 K, a c = 0.21, and 2A/T C = 4.2 ± 0.2 as so- 
lutions. We solve numerically Eqs. (0)-(|2[) to obtain the 
critical temperature T c and the zero-temperature "Mat- 
subara" gap Ag = limT^o A(itu n= Q). 

The physical gap A(T) can be obtained from A(iuj n ) 
via the analytical continuation on the real-axis |25| and 
through the relation 



A, and therefore to 2A/T C . The physical gap can be quite 
small, leading to weak-intermediate coupling phenomenol- 
ogy, even for large values of A when the electron-phonon 
coupling is balanced by a strong Coulomb repulsion, as is 
the case for fullerene compounds [Q . 

Finally, the carbon isotope coefficient ac, defined by 



Mc AT C 
AM G T c 



(7) 



is numerically evaluated by solving the ME equations (||) 
and (||) for two values of the carbon mass: Mc and Mc + 
AMc- For a single component material, it can be shown 
that isotope substitution enters the ME equations only 
through a scaling factor in the f requen cy dependence of 
the Eliashberg function, namely |2l],£6| : 



(8) 



a 2 F{Q) = T VM(2 



A(T)=Re [A(lu = A(T),T)]. 



(0) 



where M is the element mass. In a two-component sys- 
tem, like Rb 3 C6o, the electron-phonon spectral function 
a 2 F(f2) contains in principle mixed phonon modes involv- 
ing C— C, Rb— Rb and C— Rb displacements. The corre- 
sponding contributions to a 2 F(f2) scale therefore in differ- 
ent ways with isotope substitution, and one should deal 
with partial isotope effects [ pll]2^| . However, inter- molec- 
ular modes involving the alkali ions appear to couple neg- 
ligibly to electrons in A3C60 systems. Evidence for this 
conclusion comes from the negligible effect on T c upon iso- 
tope substitution of the alkali ions 127) , implying that the 
electron-phonon spectral function oFF{Q) does not con- 
tain Rb-phonon modes. In addition, the dependence of 
T c on the lattice parameter is identical both by applying 
pre ssure and by chemical substitution of the alkali atoms 
P8[ . This result points out the marginal role played by 
alkali atoms on superconductivity: they mainly tune the 
lattice constant by interspacing the buckyballs molecules 
and provide charge carriers in the conduction band. But 
they do not effectively couple to electrons. Based on this 
experimental evidence, we assume the Eliashberg function 
a 2 F(f2) is determined only by carbon modes. In this case, 
using of Eq. (||) is perfectly justified. 



In order to quantify the discrepancy between A and Aq, 
in Fig. [I] the physical and Matsubara gaps are plotted as 
function of A for a Einstein phonon model with /i = 0. 
The enhancement of A with respect to Aq is essentially 
driven by the size of the gap itself. An enlargement of 
the region relevant for RbaCgo (2A/T C — 4.2) is shown in 
inset of Fig. [I]. The discrepancy is of order of 2%, from 
2A /T C ~ 4.10 to 2A/T C ~ 4.18. The total effect is in any 
case less than the experimental error of A in RbsCeo jig] - 
Numerical solutions of Eliashberg equations in imaginary- 
axis therefore provide a quite good determination even 
of the physical gap A within the experimental accuracy 
available for RbaCgo- Note that the above discussion holds 
true even in the presence of Coulomb repulsion and for 
generic a 2 F(Q). As we have pointed out, discrepancies 
between Aq and A are essentially related to the size of 



3 Numerical analysis 

We are now in the position to analyze the experimental 
situation of RbaCgo- We solve numerically the ME equa- 
tions for different shapes of a 2 F{Q) with the constraints 
given by the measured data of RbaCeo- For the Eliash- 
berg function a 2 F(fl) we use i) a single 5-peak (Einstein 
spectrum, model I), ii) a broad spectrum defined by a 
rectangular function (model II), and iii) a broad (rect- 
angular) spectrum with an additional low frequency con- 
tribution describing the coupling to the inter-molecular 
modes (model III). 
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Fig. 2. Plot of 2A/T C , i?o> fx and /i* as functions of A as 
obtained by the numerical solutions of Eqs. (jl|)-(|2|) with an 
Einstein phonon spectrum and with the conditions T c — 30 K 
and ac = 0.21. 



3.1 Model I: Einstein phonon spectrum 

As a particularly simple but representative case we first 
consider a <5-peaked spectrum describing a single Einstein 
phonon with frequency fio and electron-phonon coupling 
constant A: 



n ). 



(9) 



The use of a single (5-function a 2 F(fi) makes the ME equa- 
tions to be characterized just by the three microscopic 
parameters A, fio and fx. Once the critical temperature 
T c = 30 K and the isotope coefficient ac = 0.21 are fixed, 
this permits to obtain a one-to-one correspondence be- 
tween the value of the ratio 2A/T C and a set of A, fio, fx. 
The numerical results are shown in Fig. |^ where we plot 
also the pseudopotential fx* obtained from the calculated 
fx via Eq. (||). 

The experimental ratio 2A/T C = 4.2 ± 0.2, together 
with T c — 30 K and ac = 0.21, is obtained by an electron- 
phonon coupling A = 3.33lQ yg, an Einstein phonon fio = 
fi\ n = 350^52 K, bare and screened Coulomb repulsions re- 



spectively fx = 1.951^^4 and M* = 0- 4 7_o:ot where the er- 
ror bars result from the experimental uncertainty of 2A/T C 
p9| . We note that the obtained value of fi\ n is compati- 
ble with the low- frequency intra- molecular modes fl30| and 
that fx* ~ 0.4 — 0.5 is close to the most accurated theo- 
retical estimations [ jBlp^ . However, A = 3.33ig;7g largely 
exceeds A ~ 1 which is the estimatation most commonly 



found in literature fij]. More importantly, as we show later, 
the value we have found is too large to prevent lattice in- 
stabilities and to ensure the validity of Migdal's theorem. 
This result does not rely on the specific shape of a 2 F(fi). 
As we show below, different shapes of a 2 F(fi) lead to sim- 
ilar results in terms of A, £2\ n and fx*. 



3.2 Model II: rectangular phonon spectrum 

Let us consider now the case of a spectral function with 
a finite width. To this end, we schematize a 2 F{fi) as a 
rectangle centered at fio and having width A fio- By using 
Eq. (0), a 2 F(fi) can therefore be written as: 



a 2 F{fi) 



\6(fi - fi + Afi /2)6(fi + Afi /2 - fi) 



21n|(/2 + AQ Q /2)/{Q - AQ Q /2)\ 



(10) 



where 9 is the Heaviside step function. The quantity AQq / Qq 
parametrizes the finite width: for Z\J?o/^o - * the single 
Einstein phonon case (model I) is recovered, while the 
limit AQq/Qq — > 2 corresponds to a constant a 2 F(Q) 
ranging from to 2i? . This latter case should be con- 
sidered merely as a mathematical limit with no physical 
relevance since, as fi — > 0, a 2 F(fi) should always vanish. 
The characteristic quantities A and fi\ n are determined as 
usual from Eqs. (ph and (^). Note that fi\ n does not coin- 
cide with fio as in the Einstein model, but it is actually 
given by: 

„-i 1/2 

)2 _ SAfi s 



fiu 



fiK 



(11) 



Hence, fi\ n is always smaller than fio and the equality 
fi\ n = fio holds true only in the limit A fio/ fio — > 0. 

The effects of the finite width of the electron-phonon 
spectrum are shown in Fig. ^, where the quantities A, fi*, 
J?i n , as well as fio and fx, obtained by requiring T c = 30 K, 
a c = 0.21, 2A/T C = 4.2 ± 0.2, are plotted as function of 
the parameter A fio/ 'fio- The error bars correspond to the 
experimental uncertainty on 2A/T C . From the analysis of 
Fig. U some important remarks can be stated. First of all, 
we see that A, fx* and fi\ n have only a weak dependence 
on the phonon spectrum width: this sustains the idea that 
T c , ac, and 2A/T C depend essentially only on the "McMil- 
lan" parameters A, fx* and fi\ n regardless of the particular 
shape of a 2 F(fi). Moreover, it shows also a lower bound 
for A (A > 2), which, although lower than the Einstein 
phonon case, is still outside the range of validity of Migdal- 
Eliashberg theory. Note that values of Afiojfio > 1.9 rep- 
resent solutions with phonon spectra exceeding ~ 2300 
K, which is the largest intramolecular phonon frequency 
[[lpO[], so that they are solutions not compatible with the 
real materials. 

Further results can be deduced from Fig. |j| In fact, 
intramolecular phonons in RbsCgo extends in a range of 
energies 393 K< fi < 2266 K If we assume that all 

these phonons couple with equal strength to the electrons, 
then this situation would correspond to fio ~ 1330 K and 
Afiojfio - 1.4. However, from Fig. |, Afi /fi = 1.4 
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Fig. 3. Plot of A, fi* , J?i n (left side axes) and flo and fi (right 
side axes) as function of the broadening of the phonon spec- 
trum Aflo I flo. T c , a and 2A/T C are fixed to their experimental 
values. 



corresponds to a phonon spectrum centered at fl ~ 550 
K, which is much lower than flo ~ 1330 K. Furthermore, 
if we keep only T c — 30 K and ac — 0.21 fixed, then 
for fl = 1330 K and Af2 /f2 = 1.4 we find 2A/T C ~ 
3.7. Again, this limiting situation is incompatible with the 
experimental data of RbsCgo- 



3.3 Model III: rectangular spectrum plus 
inter-molecular modes 

In the past, the issue of determining if intermolecular 
or intramolecular phonon modes are more responsible for 
the superconductivity phenomenon in fullerides has been 
widely debated. Although an active role of alkali-Ceo phon- 
ons has been ruled out by the zero alkali isotope effect 
p7| , intermolecular buckyball modes can in principle play 
a not negligible role in the electron-phonon coupling. In- 
deed, numerical calculations indicate that the intramolec- 
ular phonons couple relatively strongly to the conduction 
electrons [0j|j3|Q,^5|j3^] , but a contribution from very 
low-frequency intermolecular modes (librations) has been 
claimed to provide a better fit to some experimental data 
||. In literature, there are examples of numerical calcula- 
tions which favour H?! or disfavour |n||38||35| an important 
contribution of the intermolecular modes to the total A. 

In this section we consider the effects of an eventual 
intermolecular contribution to the total electron-phonon 
coupling. We schematize this contribution by adding a 
low-frequency part to the rectangular model of a 2 F(fl) 
studied above. We do not distinguish between libration (of 



typical energy 50 K g), and acoustic Cgo-Cgo modes (of 
frequency up to ~ 90— 100 K f39fl), and treat all the inter- 
molecular couplings as a continuum ranging from fl = 
to fl = Winter, where we set Anter = 100 K as the max- 
imum intermolecular frequency. The overall shape of this 
low-frequency part is of secondary importance, the only 
condition is that it should vanish as fl — > [this ensures 
that the intermolecular contribution of A is finite, see Eq. 
(H)]. We have therefore chosen the following expression for 
the total a 2 F{fl) : 

2 I Winter fl < fl < Winter 

(12) 

where, from Eq. (pj), j4i n t cr and Antra are related respec- 
tively to the intermolecular and intramolecular electron- 
phonon coupling constants, Ai n t er and Ai ntra , as: 



-4;, 



2 f2\ r . 



(13) 



and 



At 



21n|(/? + An /2)/(f2 - AQ Q /2)[ 
Moreover, form Eq. (0), the logarithmic frequency is: 



(14) 



i?l n = 



inter - 
Winter \ ^ 



nl - 



_ -^intra 

~2~X~ 



(15) 



where A = Ai ntC r + Ai ntra is the total electron-phonon cou- 
pling. In this model, the coupling to the intermolecular 
phonons is modulated just by Ai n t r, while J?i n tcr is kept 
fixed at 100 K. Of course, for Ai„t C r = the present model 
reduces to the previous model II were only intramolecular 
modes are taken into account. 

Let us start by considering the case in which all the 
intramolecular phonons (393 K< fl < 2266 K, which cor- 
responds to flo = 1330 K and Aflo/ flo — 1.4) couple with 
the same weight to the electrons. We have seen before that 
for this limiting case there is no solution compatible with 
the experimental data of RbaCgo- By switching on the in- 
termolecular electron-phonon coupling, however, it is now 
possible to find a solution which reproduces all the ex- 
perimental values of T c , a c and 2A/T C . The behaviour of 
2A/T C , together with fl\ n , /i, /i*, and the total electron- 
phonon coupling A, for fixed T c — 30 K and ac = 0.21, 
as function of Ai ntC r is plotted in Fig. [I| For Ai n t er = 
we again find 2A/T C ~ 3.7, while 2A/T C = 4.2 ± 0.2 is 
obtained for Aint, B r = l^^Qgg which corresponds to A = 



2.85™, n = 0.90 



+0.14 
-0.10' 



fi* = 0.37i°:o2 and An = 178 



+99 
-50 



K. Although now a solution exists for Ai ntC r ^ 0, the cor- 
responding values of A, fl\ n , and \i* are of the same order 
of those extracted by the previous model I and model II. 

Given Ai n t cr 7^ 0, we have also studied the effect of the 
broadening of the intramolecular modes by using differ- 
ent values of Afl§. For a given Aflo/flo, the frequency 
flo is adjusted as function of Ainter to reproduce T c = 30 
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Fig. 4. Dependence of A, [i, fj* , f2\ n and 2A/T C on the 
intermolecular electron-phonon coupling Ai nt0 r keeping fixed 
T c = 30 K and a = 0.21. 
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Fig. 5. Dependence of A, fi, n* , J?i n and J?o on the intermolec- 
ular electron-phonon coupling Ai ntor . All the data have been 
obtained by requiring T c — 30 K, ac = 0.21, and 2A/T C — 
4.2 ± 0.2 with the exception of the results for Ai nt0 r = 1.1 
where the lowest value of 2A/T C has been set equal to 4.04 
(marked by a cross) 



K, a c = 0.21, and 2A/T C = 4.2 ± 0.2. The results for 
Af2o/f2o = 1 are shown in Fig. |[ The smaller range of in- 
tramolecular phonons leads to solutions for all intermolec- 
ular couplings up to Ai n ter = 0.9. For Ai n t er = 1.1 the error 
bars of Qq becomes so large that the maximum phonon 
frequency J7 max = f?o + AQ /2 exceeds the highest possi- 
ble intramolecular phonon energy (2266 K). In particular, 
for 2A/T C = 4.04 we have obtained J? = 1930 K (marked 
by a cross in Fig. ||) which corresponds to J? max = 2895 
K. The value of J?o for 2A/T C — 4.0 would be even higher 
but too computing demanding to evaluate exactly. 

For Af2o/f2o = (Einstein phonon plus intermolecu- 
lar contribution) there are solutions also at higher values 
of Ainter, but for Aj n t r < 0.9 the resulting A, /i*, and fl\ n 
nearly exactly overlap with those of Fig. ||. Hence, com- 
pared to the case in which all the intramolecular phonons 
partecipate to the coupling [ADq/Qq = 1.4 and J?o = 
1330 K, for which we have found A = 2-85±g;||), lower 
values of A[2 /Q tend to give higher values of A. 



4 Breakdown of the adiabatic hypothesis 

We summarize in Table [j] the main results obtained by 
solving the ME equations ([!]) and (Q) under the different 
models of the el ectron-p honon spectral function a 2 F(Q) 
defined in Sees. 3T-3^3|. For model II and model III we 
report the results which give only the lower values of A 



(AQ /f2 = 1.9 for model II and Af2 /f2 = 1.4, J? = 
1330 K, and A int cr = l-47±g | 3 for model III). 

The main result which can be extracted from Table [l] 
is that, independently of the particular form of a 2 F(Q) 
considered, the experimental data of RbsCgo can be so- 
lutions of the ME equations only for very large values of 
A. Considering that model I and the lower limits of model 
II and model III reported in Table [l] correspond to quite 



Table 1. Summary of the numerical solutions of the ME equa- 
tions for T c = 30 K, a c = 0.21 and 2A/T C = 4.2 ± 0.2. For 
Model II and Model III we report only the set of values which 
defines lower limits of A. 





Oln[K\ 




A 


Model I 


350l£ 


o.47t° :£ 


q qq + 0.90 

'J-'J'J-0 .79 


Model II 


< 398ll?| 


> 0.43i°;°| 


> 2 c;i+ 81 


Model III 


> mtf 


> 0.37±° ;£ 


> 2.85±»; 8 6 3 5 
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Fig. 6. Adiabatic ratio fiph/Ep (lower line) and Migdal's pa- 
rameter P = AJ? p h/-EF (upper line) as function of the broad- 
ening of the phonon spectrum in model II. 



extreme situations, we estimate A ~ 3 ± 0.8 as the most 
plausible value for more realistic shapes of a 2 F(Q). 

Let us examine now the consequences of this result 
and the consistency with the whole ME framework. A 
first crucial point is that such high values of the electron- 
phonon coupling exceed the maximum allowed A compat- 
ible with the lattice stability. Perturbative calculations of 
the phonon propagator show that for A — > 1 the phonon 
frequencies at small momenta are renormalized to zero 
[po|l , while recent calculations on the Holstein model pre- 
dict the breakdown of the ME expansion at A > 1.25 pl[ . 
For real materials, it is argued that A ~ 1.5 is a good 
estimate of the maximum elcctron-phonon coupling com- 
patible with lattice stability Note moreover that, even 
for A well below the criterion of lattice stability, the sys- 
tem could undergo other kinds of instabilities, like charge- 
density- waves, which would prevent superconductivity. 

In addition to the problem concerning the lattice sta- 
bility, the results obtained in the last section, when an- 
alyzed together with the band structure of RbaCgo, lead 
to a serious inconsistency with respect to the validity of 
Migdal's theorem The use of Equations (|l|) and (||) 
implicitly assumes the adiabatic hypothesis according to 
which E-p 3> f2 v h, where i? p h is the typical phonon fre- 
quency and E-p is the Fermi energy. In general, the adia- 
batic hypothesis ensures that the electron-phonon vertex 
corrections, absent in the ME equation, can be neglected 
||. In fact, according to Migdal, the order of magnitude 
of the vertex corrections is at least 



P = 2 / <m 



2 F{f2) 
E F 



= A- 



fl 



E F 



(16) 



where Q ph = \ J(Wa 2 F(n) Q. Note however that the 
adiabatic condition E-p 3> fiph does not automatically 
give P < 1 since P ~ 1 could be obtained for A ^> 1. 
Hence, our results are consistent with the ME theory if 
both ilph/Ep and P are negligible. 

The electronic band structure of the fullerides is given 
by a set of very narrow subbands of width W ~ 0.5 eV 
separated from each other by gaps of order U ~ 0.5 eV 



0.6 - 

yj" 0.4 -S- 

Q. 

a 

CL 



0.2 



I i- 



4 1 



0.0 



0.0 



0.3 



0.6 



0.9 



1.2 



Fig. 7. Adiabatic ratio f2 p h/Ep (lower line) and Migdal's pa- 
rameter P = XQp^/Ep (upper line) as function of the inter- 
molecular electron-phonon coupling in model III (Af2o/f2o = 
1.0). 



or larger Q|. For the A3C60 compounds the conduction 
iiu subband is half-filled by electrons. Due to coupling 
to the phonons, the conduction electrons experience both 
intra- and inter-band scatterings. However, the inter-band 
electron-phonon couplings have only a negligible effect on 
superconductivity because they involve transitions of or- 
der W+ U 1 eV or larger JtJ . The relevant bandwidth is 
thus that of the i lu conduction band and the correspond- 
ing Fermi energy is E F = W/2 ~ 0.25 eV = 2900 K. This 
value should be compared with the characteristic phonon 
frequency scale. 

According to (|lq), in the case of an Einstein phonon 
spectrum (model I) the average phonon frequen cy J ? D h 

we 
TO. 03 

-0.02' 



reduces to H . From the previous results of Sec. ^T 
obtain therefore an adiabatic ratio fiph/Ep = 0.12 
which is moderately nonadiabatic. However, the large val- 
ue of A reported in Tabic [l] leads to an important contri- 
bution of the vertex corrections, as estimated by Eq. ([l6|), 
with a magnitude P — 0.40 j^' °2 far from to be negligible. 
The effects on P and Qp^/Ep of the broadening (model 

II) and of the inclusion of inter-molecular modes (model 

III) in a 2 F{Q) are reported in Fig. || and [?], respectively. 
For model IL the explicit expression of P can be obtained 
from Eqs. (jl(]) and ( |l6| ) and reduces to: 



P 



X 



AH 



n Q + A{2 /2\ Ep 



hi 



n - aqq/2 



(model II) (17) 



while for model III, Eq. (|l2]), P becomes: 



P 



A — Aj r 



An 



E F 



f2 + Af2 /2\ E F 



In 



n - An /2 



(model III) 



(18) 



From Fig. |^ it is apparent that the broadening affects 
only weakly the Einsteing phonon results, Af}$/ — 0, at 
least for AQq/Q^ < 1. For larger values of AQq/Qq both 
P and fiph/Ep increase signalling an even stronger viola- 
tion of the adiabatic hypothesis and of Migdal's theorem. 
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Switching on the low frequency inter-molecular couplings, 
Fig. 0, does not modify appreciably the results of Fig. 
In fact, the effect of non zero values of Ai n ter is to shift f?o 
to slightly higher frequencies than for Ainter = (see Fig. 
, leading to an almost A; n t C r-independent J7 p h- 
Summarizing the results obtained for Ep = 0.25 eV, 
we conclude that flp^/Ep > 0.12 and P > 0.4 indepen- 
dently of the shape of a 2 F(Q). Hence, the experimental 
data of RbaCgo are not consistent with the ME theory 
since both the adiabatic hypothesis, Slph/Ep <C 1, and 
Migdal's theorem, F< 1, are violated. 

How much this result depends on the precise determi- 
nation of the superconducting gap size? All over our analy- 
sis we have consider an experimental uncertainty 2A/T C — 
4.2±0.2 and we have shown that within this error bars any 
attempt to describe the experimental data breaks down 
the adiabatic hypothesis. However, recent measurements 
by surface-sensitive techniques seem to point towards a 
BCS-like value of the reduced gap 2A/T C ~ 3.53 @, and 
one could question the validity of our conclusions for such 
a small value of A. However, the violation of ME theory 
is even stronger in this case. As shown in Fig. [2] for an 
Einstein phonon spectrum (model I), BCS-like values of 
A need coupling with extremely high phonon energies in 
order to reproduce T c = 30 K and ac = 0.21. In partic- 
ular, 2A/T C — 3.6 implies J?o = 1350 K, corresponding 
to an adiabatic ratio f2 p ^/Ep ~ 0.46 and to a Migdal's 
parameter P ~ 0.49. The failure of ME theory in even 
more drastic if broader spectra (models II and III) are 
considered. Indeed, as discussed in Sec. [| in these cases a 
lowest value 2A/T C ~ 3.7 is obtained and a BCS-like gap 
2A/T C = 3.53, together with T c = 30 K and a c = 0.21, is 
incompatible with the ME framework. 



5 Discussion and conclusions 

The critical analysis of the experimental data carried on 
in the previous sections has pointed out the inconsistency 
of the standard ME theory for RbsCeo- We address now 
the origin of such inconsistency in the hypothesis that 
phonons are still the mediators of superconductivity in 
RbaCgo- To this end, let us consider some existing theoret- 
ical calculations of the coupling of the t\ u electrons to the 
H g intra-molccular phonon modes. Each phonon mode has 
frequency J?j, i = 1, . . . , 8, and couples to the ii u electrons 
via Vi = Xi/Nq, Nq being the electron density of states per 
spin at the Fermi level. In Fig. || we show a collection of 
data taken from various calculations [FP, |33| , ^ , ^5U3r| , and 
for each set of data we have estimated the corresponding 
adiabatic parameter fi^/Ep where, from Eq. (|l6|), 



(19) 



and V = 52 - Vi- The data refer to different schemes includ- 
ing tight-binding, LDA, ab-initio etc., and the error bars 
stem from the uncertainty in the calculated t\ n bandwidth 
(we have assumed E-p = 0.25 ± 0.05 eV). 
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Fig. 8. Adiabatic parameter fip^/E^ extracted from various 
calculations of the intramolecular electron-phonon interaction 
in fullerenes. 



The important point of Fig. || is that, despite of the 
large spread of the values of A/A , all these calculations 
agree in estimating the adiabatic parameter J7 p h/ 'E-p to 
be larger than about 0.4, so that Migdal's theorem breaks 
down and the whole ME framework is invalidated. Of 
course, the reason for such high values of Q-p^/Ep stems 
from the fact that, independently of details, the H g phon- 
ons have energies ranging from 30 meV to 200 meV while 
the conduction ti n electron band has a width of only 
W = 0.4 - 0.6 eV. 

With the exception of few cases (see for example Ref. 
p5|) the problem concerning the non- validity of the adi- 
abatic hypothesis is ignored or, at best, underestimated 
and the numerical results are often claimed to provide a 
strong evidence that, after all, the fullcrcne compounds are 
standard ME superconductors [0 . Our opinion is instead 
that such numerical calculations show quite clearly an ev- 
ident nonadiabaticity of the electron-phonon interaction 
for which the ME framework is completely inadeguate. 
Hence, a correct description of the superconducting state 
in fullerides should be formulated in such a way that f2 p h/ Ep 
can assume values larger than zero for which the electron- 
phonon vertex corrections, as well as finite-band effects, 
must be included from the start. As by-products of the 
very small value of Ep, also non-constant density of state 
effects and strong electron correlations may play a role 
in the superconducting properties of fullerides. However, 
these effects alone can hardly explain the high values of 
T c while, under favourable circumstances, the electron- 
phonon vertex corrections can effectively enhance the pair- 
ing interaction [^6| , leading to an increase of T c for values 
of A lower than those needed in the ME theory p7| . 

It would be therefore interesting to test whether and 
for which parameter values a theory generalized beyond 
Migdal's theorem can describe the experimental data of 
RbaCgo- The accomphshement of this task is outside the 
scope of this work. However as shown in Ref. [^3| , prelim- 
inary results suggest a positive role of the nonadiabatic 
electron-phonon effects. For example, while in the adia- 
batic ME equations a A = 1 — 4 is needed in order to 
reproduce T c = 30 K and ac = 0.21 for realistic values 
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of phonon frequencies 300 K < f2 ph < 1500 - 1800 K, a 
much more reasonable A < 1 is needed in the nonadiabatis 
theory of superconductivity where first nonadiabatic cor- 
rections are taken into account p3f . 

These results are quite interesting since they suggest 
that the experimental data of PU^Ceo could be explained 
by values of A much closer to those of the intercalated 
compounds of graphite (A ~ 0.3). This is an important 
remark because the intramolecular phonons in Ceo stem 
from the carbon-carbon bonds just as the highest phonon 
modes of graphite. In this perspective, the large discrep- 
ancy between the superconducting critical temperature in 
graphite intercalated compounds (T c ~ 0.2 K) and in ful- 
lerides would come from the difference of the electronic 
structures. Graphite compounds have large Fermi ener- 
gies locating these materials in the adiabatic regime where 
ME equations actually hold, while the small values of E-p 
of the fullcridcs give rise to the opening of nonadiabatic 
channels enhancing the pairing. 
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